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ABSTRACT 

Context. The influence of stellar activity on the fundamental properties of stars around and below 1 Mq is not well 
understood. Accurate mass, radius, and abundance determinations from solar-type binaries exhibiting various levels of 
activity are needed for a better insight into the structure and evolution of these stars. 

Aims. We aim to determine absolute dimensions and abundances for the solar-type detached eclipsing binary V636 Cen, 
and to perform a detailed comparison with results from recent stellar evolutionary models. 

Methods, uvby light curves and uvbyj3 standard photometry were obtained with the Stromgren Automatic Telescope, 
radial velocity observations with the CORAVEL spectrometer, and high-resolution spectra with the FEROS spectro- 
graph, all at ESO, La Silla. State-of-the-art methods were applied for the photometric and spectroscopic analyses. 
Results. Masses and radii that are precise to 0.5% have been established for the components of V636 Cen. The 0.85 
Mq secondary component is moderately active with starspots and Ca ll H and K emission, and the 1.05 Mq primary 
shows signs of activity as well, but at a much lower level. We derive a [Fe/H] abundance of —0.20 ± 0.08 and simi- 
lar abundances for Si, Ca, Ti, V, Cr, Co, and Ni. Corresponding solar-scaled stellar models are unable to reproduce 
V636Cen, especially its secondary component, which is ~ 10% larger and ~ 400 K cooler than predicted. Models 
adopting significantly lower mixing-length parameters I /Hp remove these discrepancies, seen also for other solar-type 
binary components. For the observed [Fe/H], Claret models for l/Hp — 1.4 (primary) and 1.0 (secondary) reproduce 
the components of V636Cen at a common age of 1.35 Gyr. The orbit is eccentric (e = 0.135 ± 0.001), and apsidal 
motion with a 40% relativistic contribution has been detected. The period is (7 = 5 270 ± 335 yr, and the inferred 
mean central density concentration coefficient, log(fc2) = —1.61 ± 0.05, agrees marginally with model predictions. The 
measured rotational velocities, 13.0 ± 0.2 (primary) and 11.2 ± 0.5 (secondary) kms~^, are in remarkable agreement 
with the theoretically predicted pseudo-synchronous velocities, but are about 15% lower than the periastron values. 
Conclusions. V636 Cen and 10 other well-studied inactive and active solar-type binaries suggest that chromospheric 
activity, and its efi'ect on envelope convection, is likely to cause radius and temperature discrepancies, which can be 
removed by adjusting the model mixing length parameters downwards. Noting this, the sample may also lend support 
to theoretical 2D radiation hydrodynamics studies, which predict a slight decrease of the mixing length parameter 
with increasing temperature/mass for inactive main sequence stars. More binaries are, however, needed for a descrip- 
tion/cahbration in terms of physical parameters and level of activity. 

Key words. Stars: evolution - Stars: fundamental parameters - Stars: activity - Stars: binaries: eclipsing - Techniques: 
spectroscopic 



1. Introduction 

As pointed o ut by e.g. Hoxie p973p, P opper ( jl997) , 
Clausen et al. (fT999l) . Torres & Ribas ^U02), Ribas ({2003)l . 
Dawson & De Robertis ([2004), and Lopez-Morales & Ribas 
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(|2005p . current stellar models, scaled to the Sun, are un- 
able to reproduce the measured temperatures and radii 
for many binary components less massive the Sun. Such 
models tend to predict too high temperatures, and radii, 
which are up to 10% too small, compared to a measured 
accuracy of about 1% for several systems. The model lu- 
minosities seem, however, to agree fairly well with obser- 
vations. Chromospheric activity has been suggested as a 
likely cause, a nd e g . Gabriel (fT969l) , Cox et al. ((TMI)) . 
Clausen et al. (|1999p . and Torres et al. (|2006p have demon- 
strated that the model fits can be improved by using a 
reduced mixing-length parameter when calculating enve- 
lope convection. The latter authors also show that there is 
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a link between activity and increased radii. Morales et al. 
(|2008p find strong radius and temperature dependencies on 
stellar activity for both single and binary low-mass stars, 
and Lopez-Morales (|2007p finds correlation between activ- 
ity level (X-ray emission level) and radius for faster rotating 
binary components, but not for slowly rotating single stars. 
On the theoretical side, Chabrier et al. p007p have, via a 
phenomenological approach, examined the consequences on 
the evolution of low-mass stars and brown dwarf eclipsing 
binaries of a) inhibiting convection due to rotation and/or 
internal magnetic field, and b) the presence of surface mag- 
netic spots. 

The current situation is summarized by Ribas et al. 
(|2008p . Clearly, accurate data for additional binaries are 
needed, as well as more models, which include first of all 
dynamo magnetic fields and the evolution of stellar rotation 
and activity (e.g. D'Antona et al. 2000). We are presently 
studying several new systems with solar-type components, 
either constant or exhibiting various levels of chromospheric 
activity (e.g. Clausen et al. I200ip . and in this paper, we 
present results from a complete analysis of the F8/G0 \{3 
system V636Cen = HD124784 = HIP69781. 

2. V636 Cen 

V636 Cen was discovered by Hoffmeister f 1958) to be a 
4^3 period eclipsing binary, and Popper (1966) reported 
sharp but single lines. Today V636Cen is known to be 
double-lined, but the lines of the secondary component are 
much fainter than those of the primary. V636 Cen is well 
detached with two quite different components in an ec- 
centric orbit (e — 0.135). Apsidal motion with a period 
of ?7 = 5 270 ± 335 yr has been detected (Clausen et al. 
I2008al hereafter CVG08). Signs of chromospheric activity 
(spots) are clearly present in the light curves, but the level 
of the light variations is modest (CVG08). Furthermore, 
we see emission from both components in the Ca ii H and 
K lines, whereas emission is not noticed; see Sect. [31 
Chromospheric emission in the Ca ii H and K lines was 
also observed by Henry et al. p996p . and X-ray emission 
was detected in the ROSAT All-Sky Survey (Voges et al. 
fT999ll . 

Our study is based on the first modern light curves 
and spectroscopic observations for V636 Cen; preliminary 
dimensions based on part of the data were derived by 
Larsen (1998) . For the spectroscopic and photometric anal- 
yses of V636 Cen presented below, we have adopted the lin- 
ear ephemeris by CVG08 (Eq. 4). Throughout the paper, 
the component eclipsed at the deeper eclipse at phase 0.0 
is referred to as the primary (p), and the other as the sec- 
ondary (s) component. 

3. Photometric elements 

The uvby light curves of V636 Cen (CVG08) contain 853 
observations in each band and were observed on 76 nights 
during six periods between March 1985 and April 1991. The 
number of observations per season are 8 (1985), 165 (1987), 
261 (1988), 233 (1989), 17 (1990), and 169 (1991). Further 
eclipse observations, but none outside eclipses, were done on 
ten nights between January 2002 and July 2007, primarily 

^ Houk p978p : spectral type of primary component. 
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Fig. 1. Ca II H and K emission in V636Cen. The vertical 
lines show the positions for the primary (full) and secondary 
(dashed) components. The FEROS spectrum of V636 Cen 
was observed at JD2451211.85 (orbital phase 0.63). 



to determine the apsidal motion period of the eccentric or- 
bit. Primary eclipse is much deeper than secondary eclipse, 
which is total and occurs near phase 0.52. 

The average observational accuracy per data point is 
about 5 mmag {vby) and 7 mmag (m), and the magnitude 
differences between the two comparison stars were found to 
be constant at that level. Additional scatter due to chro- 
mospheric activity (spots) is, however, clearly seen in the 
light curves, with total amplitudes increasing from about 
0.02 mag in y to about 0.04 mag in u, and changing some- 
what from year to year. We refer to Fig. 6 and Table 9 in 
GVG08 for further details. 

As mentioned in Sect. [H Ca ii H and K emission is 
present from both components. This is seen in two orders 
of two FEROS spectra from February 1999 (Table [6]), and 
(dark) spots may therefore be present at the surface of both 
stars. Fig.[T]illustrates that the emission features of the two 
components appear to have about the same strength, but 
since the primary dominates the total light of V636 Cen, the 
real emission from the secondary is by far the strongest. 

Between 1987 and 1991, the depths of the total sec- 
ondary eclipses were nearly identical, whereas those of the 
primary eclipses differed significantly. During this period, 
the activity must therefore have been dominated by con- 
tributions from the secondary component. As mentioned 
below, this is supported by frequency analyses of the out-of- 
eclipse data. On the other hand, the depths of both eclipses 
varied from 2002 to 2006. This is not due to the slow ap- 
sidal motion detected for V636Cen, and must be caused 
by spots at both stars during at least some of those years. 
The RossbjO numbers for the primary and secondary com- 
ponents are approximately 0.44 and 0.10, respectively. This 
places the secondary in the range where other stars tend to 
show spot activity and the primary close to onset of spot- 
tedness (Hall lT99t| . 

^ Defined as the ratio of the rotation period to the convective 
turnover time. 
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Table 1. Photometric solutions for V636Cen from 
JKTEBOP analyses of the 1987 and 1988 hght curve ob- 
servations (426 points per band). 



y h V u 



i (°) 


89 65 


89 62 


89 62 


89 63 


±7 


±7 


±8 


±9 


' P 


0746 


0748 


0750 


0751 




±5 


±5 


±6 


±6 




0613 


0614 


0613 


0614 




±3 


±3 


±3 


±3 


k 


0.822 


0.822 


0.818 


0.817 




±5 


±5 


±6 


±6 


Tp + Ts 


0.1359 


0.1362 


0.1363 


0.1365 




±6 


±5 


±7 


±7 


Up 


0.59 


0.68 


0.77 


0.81 


Us 


0.73 


0.84 


0.95 


0.99 


Vp 


0.37 


0.43 


0.49 


0.58 


Vs 


0.44 


0.51 


0.59 


0.70 


fJ s / J p 


0.417 


0.370 


0.286 


0.236 


±10 


±9 


±8 


±7 


Ls 1 Lp 


0.265 


0.233 


0.176 


0.144 


±3 


±3 


±3 


±3 


a (mmag.) 


11.3 


11.4 


13.2 


14.5 



Note 1: Orbital parameters of e = 0.1348 and co = 281?60 were 
assumed (Table[5ll, and linear limb darkening coefficients by Van 
Hamme ([1993P were adopted. 

Note 2: The errors quoted for the free parameters are realistic 
values determined from 1 000 Monte Carlo simulations, e and 
u were perturbed corresponding to their uncertainties, and the 
limb darkening coefficients within ±0.05. 

Table 2. WD solutions from the 1987 light curve obser- 
vations (165 per band). 





V 


h 


V 


u 




89.80 


89.63 


89.83 


89.81 


±3 


±2 


±2 


±4 




14.504 


14.465 


14.395 


14.397 




±22 


±33 


±38 


±49 


ns 


14.738 


14.666 


14.725 


14.720 




±17 


±34 


±40 


±56 


rp' 


0.0737 


0.0739 


0.0743 


0.0743 




0.0601 


0.0604 


0.0601 


0.0602 


k 


0.815 


0.817 


0.809 


0.810 


Tp + rs 


0.1338 


0.1343 


0.1344 


0.1345 


Ts (K) 


4965 


4984 


4916 


4962 


±6 


±5 


±7 


±7 


Ls / Lp 


0.286 


0.256 


0.199 


0.172 


Lo, n 


12.2 


15.0 


17.2 


31.8 


±3.7 


±3.2 


±3.3 


±6.6 


ARi (°) 


23.9 


25.7 


29.6 


31.8 


±1.0 


±0.8 


±1.0 


±1.2 


L02 n 


233.0 


235.3 


233.8 


255.6 


±7.6 


±3.7 


±3.5 


±10.1 


AR2 C) 


17.7 


19.3 


22.5 


25.7 


±1.4 


±0.8 


±1.2 


±1.2 


a (mmag.) 


6.2 


6.2 


7.7 


9.1 



^ relative volume radius 



Note 1: Tp was assumed to be 5900 K; for orbital parameters 
and limb darkening coefEcients see Table [1] 
Note 2: The errors quoted for the free parameters are the 
formal probable errors determined from the iterative least 
squares solution procedure. 



Table 3. WD solutions from the 1988 light curve obser- 
vations (261 per band). See Tabled) 





y 


I) 






i C) 


89.56 


89.55 


89.56 


89.40 


±3 


±3 


±3 


±2 




14.495 


14.470 


14.455 


14.477 




ztzo 




ztZO 






14.656 


14.639 


14.707 


14.623 




4-07 


±zO 


4-1/1 
±14 


±iy 




0.0738 


0.0739 


0.0740 


0.0739 




0.0604 


0.0605 


0.0602 


0.0606 


k 


0.818 


0.819 


0.814 


0.820 


Tp ± Ts 


0.1342 


0.1344 


0.1342 


0.1345 


Ts (K) 


4940 


4959 


4897 


4914 


±4 


±4 


±4 


±6 


T IT 

^s 1 l^p 




U.zoU 


U.iyo 


U.loo 


Lo, n 


107.9 


118.8 


116.5 


122.0 


±1.6 


±1.2 


±1.4 


±2.6 


ARr (°) 


25.9 


28.7 


31.8 


31.1 


±0.4 


±0.4 


±0.4 


±0.6 


LO2 (°) 


202.4 


212.1 


209.9 


214.9 


±1.7 


±1.6 


±1.6 


±2.4 


AR2 (°) 


25.7 


23.8 


28.6 


29.9 




±0.6 


±0.7 


±0.7 


±0.9 


a (mmag.) 


6.5 


5.9 


6.8 


9.7 



Out-of-eclipse phases were covered on several nights 
during 1987, 1988, 1989, and 1991. Frequency analyses of 
the corresponding data, done separately for each year and 
band, reveal the existence of one significant period per year. 
The results from the four bands agree well, and the average 
periods and estimated uncertainties are 4'?05±0^06 (1987), 
4'?16 ± Qt^Ql (1988), 4"?13 ± 0^01 (1989), and 4^06 ± 0^12 
(1991). They are all shorter than the orbital period 4'?28, 
but transformed to equatorial surface velocities they agree 
quite well with the observed rotational velocities of both 
components (Sect.S]). 

Based on this information on eclipse depths and period- 
icity of the out-of-eclipse variability, we conclude that the 
1987-1991 surface activities can be modeled adequately by 
including a few dominating spots (or spot groups) on the 
surface of the secondary component. 

We have used the Wilson-Devinney (WD) model and 
code for the analyses. As described below, some initial 
analyses, which ignore surface activity, were based on the 
JKTEBOP program. We refer to Clausen et al. (|2008b[ 
hereafter CTB08) for references and details on the binary 
models and codes, and on the general approach applied. 
In tables and text, we use the following symbols: i orbital 
inclination; r relative radius; k = rs/rp] surface poten- 
tial; u linear limb darkening coefficient; y gravity darkening 
coefficient; J central surface brightness; L luminosity; Tog 
effective temperature. 

A linear limb darkening law has been assumed through- 
out with coefficients adopted from Van Hamme (|1993p 
and/or Claret (200(7). They differ by 0.02-0.10 only, lead- 
ing to identical photometric elements within errors. Gravity 
darkening coefficients/exponents corresponding to convec- 
tive atmospheres were applied, and bolometric reflection 
albedo coefficients of 0.5 were chosen in the WD analy- 
ses, again due to convection. For JKTEBOP, the simple 
bolometric reflection model was used. The mass ratio be- 
tween the components was kept at the spectroscopic value 
[Ms/Mp = 0.812±0.002). In order to model the spots prop- 
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Fig. 2. Fit of the theoretical light curve of the WD solution from Table [2] to the 1987 y observations. 




Phase 

Fig. 3. Fit of the theoretical light curve of the WD solution from Table [3] to the 1988 y observations. 



erly, the (very nearly spherical) components were assumed 
to rotate at the periods determined from the frequency 
analyses mentioned above. As initial analyses showed that 
the orbital eccentricity is not as well constrained from the 
light curve analyses as from the spectroscopic orbit, we have 
adopted the spectroscopic result throughout; see Table [S] 
The longitude of periastron uj is then defined through 
e cos UJ from the well determined phase of central secondary 
eclipse, which did not change significantly relative to pri- 
mary eclipse during the light curve observations from 1987 
to 1991. 

As the activity level of V636 Cen is modest, it is ap- 
propriate to first analyse the light curves without tak- 
ing spots into account, and to adopt a simple but ad- 
equate binary model/code like JKTEBOP, which allows 
easy determination of realistic errors of the photometric 
elements from Monte Carlo simulations. Since the average 
out-of-eclipse level of the 1991 observations is significantly 
off, and since there are very few 1985 and 1990 observa- 
tions, these data sets were not included in the analyses. 
Results from the combined 1987 and 1988 observations, 
which both cover all phases well, are given in Table [1] 
Including also the 1989 observations increases the scatter 



significantly but leads to identical solutions within errors. 
As seen, the results from the four bands agree remarkably 
well, and relative radii with a realistic precision of about 
1% are obtained. Changing the adopted orbital eccentricity 
by -1-0.005, i.e. four times the error of the spectroscopic re- 
sult, decreases the relative radii of both components by only 
0.5%. Adjusting e and uj leads to spurious results, due to 
asymmetries in the observed light curves caused by spots. 
As expected, photometric elements from similar WD anal- 
yses, not accounting for spots, agree very well with those 
from JKTEBOP. 

Results from WD analyses of the 1987 and 1988 obser- 
vations are presented in Tables [2] and [3l and out-of-eclipse 
comparisons in the y band are shown in Figs. [2] and [3l The 
activity is modeled by including two rather big spots at 
the surface of the secondary component. Most probably, 
this simplified approach does not represent the real surface 
structure in detail, but it reproduces the light variations 
quite closely, and somewhat better than if only one spot is 
included. Based on extensive tests, the spots were placed at 
the equator of the secondary component, and a ratio of 0.9 
between the temperatures of spotted and unspotted areas 
was selected. Longitudes {LO) and angular radii {AR) of 
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Table 4. Adopted photometric elements for V636 Cen. 



89? 65 ±0? 10 
0.0740 ± 0.0003 
0.0603 ± 0.0003 



y 

0.269 

±15 



0.237 

±15 



0.178 

±15 



0.152 

±15 



Note: The individual luminosity ratios are based on the mean 
stellar and orbital parameters and include the average activity 
level (spots) in 1987 and 1988. 



the spots were adjusted together with orbital and stellar 
parameters. 

As seen, the eight independent WD solutions yield 
nearly identical relative radii of the components, which are 
systematically about 1.5% smaller than the JKTEBOP re- 
sults (Tabled]). The orbital inclinations of the 1987 solu- 
tions are on average about 0?25 higher than those from the 
1988 solutions, but tests show that the realistic error of i 
is at least 3-4 times larger than the formal errors listed. 
The luminosity ratios given in Tables [2] and [Sj which for 
the sake of comparison do not include the spots, agree well 
within realistic errors; we will return to this point below. 
For each year, the spot parameters from the four bands dif- 
fer slightly, but not more than can be expected from the 
simple approach we have adopted and the limited number 
of observations available. The two spots applied cover about 
10% of the surface of the secondary component. Finally, the 
rms errors (a) of the fits to the observations are close to 
the observational errors per data point. 

Similar WD analyses of the 1989 observations were not 
successful. The ascending branch of secondary eclipse is not 
covered, and the level of activity seems to have changed dur- 
ing the last half of the observing period. Secondary eclipse 
was not covered in 1991, and very few 1985 and 1990 ob- 
servations exist, so these data sets were not analysed. 

As a further test of the simple spot model adopted, it 
is of interest to compare the 1987 and 1988 light curves 
cleaned for spot effects. This is easily done by calculating, 
year by year and band per band, theoretical light curves 
with and without spots and then subtract the difference be- 
tween them from the observations. First, the cleaned 1987 
and 1988 light curves agree well, and second, JKTEBOP 
analyses of the combined cleaned data yield stellar radii, 
which are now in perfect agreement with those from the 
WD solutions. The orbital inclination is close to the mean 
of the 1987 and 1988 WD resuhs. 

The adopted photometric elements are presented in 
Table The main result of the photometric analyses is 
that we have, in spite of the surface activity, been able to 
determine very accurate relative radii for the components 
of V636 Cen. The system is well detached and the stars are 
quite different, especially with respect to temperature and 
luminosity. As seen, the orbital inclination is high, and con- 
sequently secondary eclipse is total. The luminosity ratios 
listed in Table |4] are based on the mean stellar and orbital 
parameters and include the average spot level seen in the 
1987 and 1988 observations; errors reflect the difference be- 
tween the two years. This approach was adopted in order 
to allow calculation of realistic uvby indices for the compo- 
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Fig. 4. Spectroscopic orbital solution for V636 Cen (solid 
line: primary; dashed line: secondary) and radial velocities 
(filled circles: primary; open circles: secondary). The dot- 
ted line (upper panel) represents the center-of-mass veloc- 
ity of the system. Phase 0.0 corresponds to central primary 
eclipse. 



nents; see Table |S1 Due to the spots on the surface of the 
secondary component, the ratio of the light contributions 
outside eclipses of course change significantly through the 
orbit. However, the b — y, mi, and ci indices of the compo- 
nents remain nearly unaffected. 



4. Spectroscopic orbit 

The radial velocity measurements of V636 Cen were ob- 
tained with the photoelectric cross-correlation spectrom- 
eter CORAVEL (Baranne et al. [T9791 Mayor 1985) oper- 
ated at the Danish 1.5-m telescope at ESO, La Silla. In 
total, 38 and 39 measurement were obtained for the pri- 
mary and secondary components, respectively, on 30 nights 
between March 1986 and June 1989. The internal veloc- 
ity errors range between 0.39 and 0.65 kms~^ (median 
0.44 kms~^) for the primary component and between 0.71 
and 1.81 kms"""^ (median 0.88 kms~^) for the fainter sec- 
ondary component. The observations, which cover the out- 
of-eclipse phases well, are listed in Table lA.ll 

The spectroscopic orbits were calculated using the 
method of Lehman-Filhes implemented in the SBOF0 pro- 
gram (Etzel 1^041) . which is a modified and expanded ver- 
sion of an earlier code by Wolfe, Horak & Storer (1967|. 
Following several tests, we have decided to weight the in- 



^ Spectroscopic Binary Orbit Program, 
http://mintaka.sdsu.edu/faculty/etzei7] 
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Table 5. Spectroscopic orbital elements of V636Cen. 





Pri TTifirv 

J. ± illlcLl. y 




PviTTl FlTV-l-SipPnTl H aT'V 
A ± iiiicti. y 1 kj t^v., W11LJ.CX1. y 

Adopted solution 


(kms^^) 


73.54 ± 0.07 




73.53 ±0.10 


Ks (kms-^) 




90.48 ± 0.25 


90.51 ±0.19 


7 (kms~"^) 


-14.23 ±0.07 


-14.31 ± 0.24 


-14.25 ±0.07 


e 


0.1353 ±0.0010 


0.1339 ±0.0029 


0.1348 ±0.0012 


^ n 


281.57 ±0.58 


281.23 ± 1.78 


281.60 ±0.63 


T-2 400 000=^ 


47223.1225 ± 0.0066 


47223.1200 ± 0.0199 


47223.1233 ± 0.0072 


ap^ (kms-i) 


0.33 




0.42 


as"" (kms-^) 




1.23 


0.79 



Time of periastron 
^ Standard deviation of a single radial velocity 

Note: An orbital period of 4"?2839423 has been assumed (CVG08). 



dividual radial velocities according to the inverse square 
of their internal errors. The SBl and SB2 solutions for 
V636 Cen are given in Table [5l and the observations and 
computed orbit are shown in Fig. |4l As seen, the SB2 solu- 
tion, which we adopt, agrees well with the individual SBl 
solutions. The semiamplitudes yield minimum masses ac- 
curate to about 0.4% for both components, and the eccen- 
tricity and periastron longitude of the orbit are well con- 
strained. Also, the standard deviations of a single velocity 
agree well with the internal median velocity errors. 

The measured vsmi values are 13.0 ± 0.2 kms~^ and 
11.2 ± 0.5 kms^^, respectively, with the velocity of the 
secondary being more uncertain due to its much weaker 
lines. These rotational velocities are supported by the fre- 
quency analyses of the out-of-eclipse photometry described 
in Sect 121 As discussed in Sect. ISl the rotation of the com- 
ponents corresponds closely to pseudo-synchronization. 



5. Abundances 

For abundance determinations, we have obtained two high- 
resolution spectra of V636Cen with the FEROS fiber 
echelle spectrograph at ESO, La Silla in February 1999; see 
Table IHl Details on the spectrograph and the reduction of 
the spectra are given by CTB08. The late-type components 
of V636 Cen have rich absorption line spectra, and unlike 
many known solar-type systems fTable lTUl) . the components 
of V636 Cen are relatively slow rotators. This is a clear ad- 
vantage for the abundances studies. The ideal case would 
be to have enough spectra for disentangling, since mu- 
tual blending is avoided and the reconstructed component 
spectra get a higher signal-to-noise ratio than each of the 
observed spectra. However, as will be demonstrated below, 
the wide spectral coverage of FEROS provides a large num- 
ber of unblended lines from both components which are us- 
able for the analyses of the combined spectra. Furthermore, 
the two spectra have been observed at different phases giv- 
ing opposite line shifts. Since the primary component is 3- 
5 times more luminous than the secondary, depending on 
wavelength, its lines are only diluted by a factor of about 



About 10 spectra with S/N ~ 100, well distributed in phase, 
would be required 



Table 6. Log of the FEROS observations of V636 Cen. 



HJD-2 400 000^ 


phase 


t ' 




51210.8718 


0.4017 


2100 


160 


51211.8484 


0.6296 


2400 


200 



" Refers to mid-exposure 
' Exposure time in seconds 
Signal-to-noise ratio measured around 6070 A 



1.3, meaning that an accuracy close to that for a similar sin- 
gle star with spectra of comparable quality can be reached. 
For the secondary, the situation is on the other hand worse; 
its dilution factor is about 4.7. A lower accuracy must there- 
fore be expected for the secondary due to higher errors of 
the intrinsic equivalent widths. In general, also intrinsically 
stronger lines, for which the equivalent widths arc less sen- 
sitive to abundance and more sensitive to e.g. microtur- 
bulence and pressure broadening, have to be used for the 
secondary. This is in principle avoided in disentangling, but 
the signal-to-noise ratio of the reconstructed spectrum for 
the secondary would on the other hand still be much lower 
than for the primary. 

The basic approach followed in the abundance analyses 
is described by CTB08. The versatile VWA tool, now ex- 
tended to analyses of double-lined spectra, was used; we re- 
fer to Bruntt et al. ( 2004 , 2008J for an updated description 
of VWA. It uses the SYNTH software (Valenti & Piskunov 
il996|) to generate the synthetic spectra. Atmosphere mod- 
els were interpolated from the grid of modified ATLAS9 
models by Heiter et al. (|2002[) . Line information was taken 
from the Vienna Atomic Line Database (VALD; Kupka et 
al. I1999P , but in order to derive abundances relative to the 
Sun, \og{gf) values have been adjusted in such a way that 
each measured line in the Wallace et al. (|1998p solar at- 
las reproduces the atmospheric abundances by Grevesse & 
Sauval (|1998p . Analyses of a FEROS sky spectrum repro- 
duce these adjustments closely. 

The abundance results derived from all useful lines in 
both spectra are presented in Table H Except for a few 
ions, we have only included lines with measured equivalent 
widths above 10 mA and below 80 mA (primary) and 50 
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Table 7. Abundances ([El./H]) for the primary and sec- 
ondary components of V636 Cen determined from the two 
FEROS spectra. 







Primary 






Secondary 




Ion 


[El./H] 


rms 


Nt/N, 


[El./H 


] rms 


Nt/N, 


Si I 


-0.17 


0.12 


15/15 


-0.26 


0.15 


3/3 


Cal 


-0.15 


0.11 


11/9 


—0.10 


0.23 


12/9 


Til 


-0.19 


0.15 


5/5 


-0.34 


0.16 


8/6 


Tin 


-0.22 


0.10 


6/4 








Vi 


-0.27 


0.04 


6/4 


-0.15 


0.17 


13/12 


Cri 


-0.19 


0.15 


6/4 








Crn 


-0.05 


0.16 


6/5 








Fei 


-0.19 


0.09 


136/110 


-0.11 


0.19 


48/46 


Fen 


-0.20 


0.08 


16/12 








Co I 


-0.15 


0.15 


7/5 








Nil 


-0.22 


0.13 


39/25 


-0.19 


0.17 


12/9 



Note: Nt is the total number of lines used per ion, and N; is 
the number of different lines used per ion. Ions with at least 3 
lines measured are included. 



mA (secondary). Comparing the results from the individ- 
ual spectra, we find no significant differences. The effec- 
tive temperatures, surface gravities and rotational veloci- 
ties listed in Table [5] were adopted. Microturbulence veloc- 
ities were tuned until Fe i abundances were independent 
of line equivalent widths, and the resulting values are 1.15 
(primary) and 1.35 (secondary) kms~^. The calibration by 
Edvardsson et al. (flMg]) predicts 1.25 ±0.31 kms-^ for the 
primary component, whereas the secondary lies outside its 
effective temperature range. For the adopted effective tem- 
peratures we see no dependency of the abundance on exci- 
tation potential, which, however, occurs if they are changed 
by more than ±100 - 150 K. 

As seen, a robust [Fe/H] is obtained for the primary, 
with identical results from Fe I and Fe ii lines. Within er- 
rors, the less precise result for the secondary from Fe i lines 
agrees; unfortunately no suitable, unblended Fe ii lines are 
available. Changing the primary model temperatures by 
±100 K modifies [Fe/H] from the Fe i lines by ±0.05 dex, 
and a similar but opposite effect is seen for Fe ii lines. If 
0.3 kms~^ higher microturbulence velocities are adopted, 
[Fe/H] decreases by about 0.05 dex for both neutral and 
ionized lines. Finally, changing the adopted luminosity ra- 
tio between the components (Table 2]) by +1 cr changes 
[Fe/H] by +0.015 dex for the primary and -0.04 dex for the 
secondary. Taking these contributions to the uncertainties 
into account, and giving higher weight to the primary, we 
adopt [Fe/H] = -0.20 ±0.08 for V636Cen. Except for Cr ii 
we also find relative abundances close to -0.20 dex for the 
other ions listed in Table [71 including the a-elements. The 
reason for the Cr ii discrepancy is not clear, but a simi- 
lar result was recently found for the GO V system WZ Oph 
(CTB08). 

We have also done abundance analyses based on the re- 
cent grid of MARCS model atmospheres (Gustafsson et al. 
I2008p . which adopt the solar composition by Grevesse et 
al. (|2007p . Adjusting the log{gf) accordingly, as described 
above, we obtain abundances for V636Cen which agree 
with those listed in Table [7] within about ±0.05 dex. 

As a supplement to the spectroscopic abundance anal- 
yses, we have derived metal abundances from the de- 
reddened uvby indices for the individual components (see 
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Table 8. Astrophysical data for V636 Cen. 





Primary 


Secondary 


A riQrili i"i"p (1 i tn pn m nn q ■ 






M/Mq 


1.052 ± 0.005 


0.854 ± 0.003 


R/Rq 


1.018 ±0.004 


0.830 ± 0.004 


log 5 (cgs) 


4.444 ± 0.004 


4.532 ± 0.005 


V sini^ (kms""*^) 


13.0 ±0.2 


11.1 ±0.5 


Vsync^ (kmS~^) 


12.0 ±0.1 


9.8 ±0.1 


Vpsync (kmS ) 


13.3 ±0.1 


10.9 ±0.1 


Vperi'^ (kmS""*") 


15.9 ±0.1 


13.0 ±0.1 


Photometric data: 






V 


8.963 ±0.017 


10.389 ± 0.049 


ih-vY 


0.382 ± 0.004 


0.520 ± 0.009 


mi' 


0.179 ±0.007 


0.357 ±0.017 




0.315 ±0.018 


0.169 ±0.020 


E{b - y) 


0.005 ± 0.008 




5900 ± 85 


5000 ± 100 


DOi 


4.61 ±0.06 


5.78 ±0.09 


log L/Lq 


0.05 ±0.03 


-0.41 ±0.04 


BC 


-0.06 


-0.30 


Mv 


4.67 ±0.06 


6.08 ±0.09 


V -Mv 


4.29 ±0.07 


4.31 ±0.10 


Vo - Mv 


4.27 ±0.08 


4.29 ±0.11 


Distance (pc) 


71.4 ±2.5 


72.2 ±3.6 


Abundance: 






[Fe/H] 


-0.20 


±0.08 


Observed rotational velocity 



^ Equatorial velocity for synchronous rotation 

Equatorial velocity for pseudo-synchronous rotation 
Refers to periastron velocity 
° Not corrected for interstellar absorption/reddening 

Note: Bolometric corrections {BC) by Flower l|1996p have been 
assumed, together with T^^^q = 5780 K, BCq = —0.08, and 
M,„,0 = 4.74. 



Table [HI) and the calibration by Holmberg et al. (|2007p . For 
the primary component, the result, [Fe/H] = —0.23 ± 0.11, 
agrees perfectly well with the spectroscopic determination, 
whereas an unrealistically low value, [Fe/H] — — 0.77±0.17, 
is obtained for the secondary. Their previous calibration 
for late type stars (Nordstrom et al. I2004p gives [Fe/H] = 
—0.86 ± 0.20. The quoted [Fe/H] errors include the uncer- 
tainties of the photometric indices and the published spread 
of the calibrations, and the reason for the discrepancy be- 
tween the primary and secondary results is not clear. It 
may, however, very likely be due to mi deficiency for the 
secondary related to activity, as in fact seen for active bi- 
naries (Gimenez et al. 119911) . rather than a strong metal 
underabundance. Adopting the definitions by Gimenez et 
al., we find 5mi ~ 0.11 and 5ci ^ 0.00 for the secondary. 
If we modify mi accordingly and assume that 6 — ?/ is 
unaffected, the two calibrations mentioned above lead to 
[Fe/H] = -0.34 ± 0.18 and [Fe/H] = -0.36 ± 0.21, respec- 
tively, i.e. reasonable agreement with the result for the pri- 
mary is reached. 
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6. Absolute dimensions 

Absolute dimensions for the components of V636 Cen are 
calculated from the elements given in Tables 2] and \5\ As 
seen in Table HI both masses and radii have been obtained 
to an accuracy of about 0.5%. 

Individual standard uvby indices are included in 
Table [51 as calculated from the combined indices of 
V636 Cen outside eclipses (CVG08) and the luminosity ra- 
tios (Table [4]). The calibration by Olsen ,1988, the uvby 
indices of the primary component and P observations ob- 
tained during the total secondary eclipse {f3 = 2.601±0.008) 
then yield a very small interstellar reddening of E{b — y) — 
0.005 ±0.008. 

The adopted effective temperatures (5900 K, 5000 K) 
were calculated from the calibration by Holmberg et al. 
(|2007| . assuming [Fe/H] = -0.20 (Sect. [5]). The uncertain- 
ties include those of the uvby indices, E{b — y), [Fe/H] and 
the calibration itself. Temperatures based on the calibra- 
tions by Alonso et al. (|1996p and Ramirez & Melendez 
(|2005p agree within errors, but are, however, systematically 
about 100 K lower for the primary component. 

From the orbital and stellar parameters, the observed 
rotational velocities, and the apsidal motion period deter- 
mined by CVG08, we derive a mean central density coeffi- 
cient of log(fc2) = —1.61 ±0.05. The observed rate of change 
of the longitude of periastron cu = 0.00080±0.00005 °/c was 
corrected for a 40% relativistic contribution of a; = 0.00032 
°/c, as calculated from general relativity (Gimenez. I1985p . 
The higher order and terms are negligible. From the- 
oretical coefficients, calculated for the models described in 
Sect. O along the Hues described by Claret (fTOMl [TM7|) . 
we derive a theoretical log(fc2) of —1.74. 

As seen in Table |S1 the measured rotational velocities 
are in remarkable agreement with the theoretically pre- 
dicted pseudo-synchronous velocities (Hut il981i. Eq. 42), 
but are about 15% lower than the periastron values. The 
turbulent dissipation and radiative damping formalism of 
Zahn (|1977[ 1989) predicts synchronization time scales of 
0.1 Gyr for both components and a time scale for circular- 
ization of 7.5 Gyr, whereas the hydrodynamical mechanism 
by Tassoul & Tassoul (|19971 and references therein) yields 
much lower values of 0.3 Myr (primary synchronization), 
0.4 Myr (secondary synchronization), and 0.04 Gyr (circu- 
larization). We refer to Claret & Cunha (|1997p and Claret 
et al. H1995P for details on the calculations. From model 
comparisons (Sect. [7?2l), we derive an age for V636Cen of 
about 1.35 Gyr, meaning that both theories are in agree- 
ment with the observed pseudo-synchronization. On the 
other hand, the mechanism by Tassoul & Tassoul predicts 
a circular orbit, which is certainly not the case. 

The distance to V636 Cen was calculated from the 
"classical" relation (see e.g. CTB08), adopting the solar 
values and bolometric corrections given in Table [8] and 
Av/E{b — T/) = 4.27, and accounting for all error sources. 
As seen, nearly identical results are obtained from the two 
components. The mean distance, 71.8 pc, which has been 
established to about 4%, is in perfect agreement with the re- 
sult of the new Hipparcos reduction by van Leeuwen p007p . 
72.2 ± 4.7 pc, but is marginally larger than the original 
Hipparcos result 65.1 ± 4.7 pc (ESA IWf). We note that 
V636 Cen belongs to the group of eclipsing binaries within 
125 pc, discussed by Popper (|1998p . which could be useful 
for improving the radiative flux scale. 




6000 



5500 



5000 



Fig. 5. V636Cen compared to models calculated for 
the measured abundance [Fe/H] = —0.20. Tracks for the 
component masses (full drawn, thick) and isochrones for 
0.5, 2.0, 4.0, and 6.0 Gyr (dashed) are shown. The uncer- 
tainty in the location of the tracks coming from the mass 
errors are indicated (dotted lines). To illustrate the effect of 
the abundance uncertainty, tracks for [Fe/H] = —0.12 are 
included (full drawn, thin). 



7. Discussion 

In the following, we first compare the absolute dimensions 
obtained for V636 Cen with properties of recent theoretical 
stellar evolutionary models. We then compare V636 Cen to 
other well-studied eclipsing binaries with at least one com- 
ponent in the 0.8-1.1 Mq mass interval. 



7.1. Comparison with solar-scaled models 

Figs. [51 El and[71illustrate the results from comparisons with 
the Yonsei-Yale {Y'^) evolutionary tracks and isochrones 
by Demarque et al. (I2004|FI . The mixing-length parame- 
ter in convective envelopes is calibrated using the Sun, 
and is held fixed at I /Hp — 1.7432. The enrichment law 
Y = 0.23 + 2Z is adopted, together with the solar mix- 
ture by Grevesse et al. (I1996p . leading to {X,Y,Z)q = 
(0.71564,0.26624,0.01812). Only models for [a/Fe] = 0.0 
have been considered. We refer to CTB08 for a brief de- 
scription of other aspects of their up-to-date input physics. 

As seen from Fig. [5l models for the observed masses 
and abundance, [Fe/H] = —0.20, equivalent to {X,Y,Z) = 
(0.73475,0.25350,0.01175), are significantly hotter than ob- 
served, especially for the secondary component. The un- 
certainty of [Fe/H] is ±0.08 dex, and models for [Fe/H] = 
-0.12, equivalent to {X,Y,Z) = (0.72809,0.25790,0.01397), 
marginally fit the primary but not the secondary compo- 
nent. In order to fit the primary perfectly, Y^ models with 
a heavy element content close to solar, or a higher He abun- 
dance, are required. Such models would, however, still be 
unable to reproduce the secondary. The same picture is 



jhttp : //www . astro . yale . edu/demarque/yystar . html| 
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Fig. 6. V636 Cen compared to models calculated for 
[Fe/H] = -0.20. Isochrones for 0.5, 2.0, 4.0, and 6.0 Gyr 
are shown. For comparison, 2 Gyr isochrones for [Fe/H] — 
-0.12 (dashed) and [Fe/H] = 0.00 (dotted) are included. 
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Fig. 7. V636 Cen compared to models calculated for 
[Fe/H] = -0.20. Isochrones for 0.5, 2.0, 4.0, and 6.0 Gyr 
are shown. For comparison, 2 Gyr isochrones for [Fe/H] = 
-0.12 (dashed) and [Fe/H] = 0.00 (dotted) are included. 

seen for e.g. Victoria-Regina (Vand enBerg et al., [2006)FI 
and BaSTI (Pietrinferni et al.. [MMlFI models, which differ 
slightly from Y^, e.g. with respect to input physics and He 
enrichment law. We refer to CTB08 for a brief comparison 
of the three model grids. 

The most direct comparison, based on the scale inde- 
pendent masses and radii of V636 Cen, is shown in Fig. [S] 
For the primary component, the models predict an age of 
2.6-3.2 Gyr, depending on the assumed [Fe/H], whereas the 

^ [littp : //wwwl . cadc-ccda.hia-iha.nrc-cnrc .gc . ca/cvo/| 
community/VictoriaReginaModels/ 
^ |http : //www . te . astro . it/BASTI/index .php| 



Table 9. Ages (Gyr) for the primary (p) and secondary (s) 
components of V636 Cen determined from comparison be- 
tween their masses and radii and Claret models for different 
mixing length parameters and abundances. 





[Fc/H] 




l/H,{s) 


Age(p) 


Age(s) 


1 


-0.20 


1.40 


1.00 


1.33 ±0.13 


1.37 ±0.55 


2 


-0.20 


1.50 


1.10 


1.82 ± 0.13 


3.10 ±0.55 


3 


-0.12 


1.60 


1.10 


2.46 ±0.12 


3.25 ±0.53 


4 


-0.20 


1.68 


1.68 


2.40 ±0.13 





Note: Ages in Gyr. Except for those with the solar I /Hp = 
1.68 (case 4), the models fit the observed eflFective temperatures 
within errors. See text for details. 



secondary is predicted to be more than two times older. At 
the predicted age of the primary, the observed radius of 
the secondary, which has been determined to an accuracy 
of 0.5%, is about 10% larger than that of the corresponding 
model. Mass-luminosity comparisons are shown in Fig. [71 
Within the errors of the observed luminosities, the models 
predict identical - but lower - ages for the two components. 

In conclusion, the secondary component of V636 Cen is 
significantly larger and ~ 400 K cooler than predicted by 
models which adopt a mixing-length parameter matching 
the Sun. For the observed [Fe/H], such models are ^ 200 K 
hotter than the primary component, indicating a discrep- 
ancy as well. 

7.2. Comparison with mixing-length "tuned" models 

As mentioned in Sect.[Tl several authors have demonstrated 
that models which adopt a reduced envelope convection 
mixing-length parameter fit active late-type binary com- 
ponents better. Since V636 Cen exhibits intrinsic variation 
and Ca II H and K emission (Fig. [T]) , and has relatively high 
rotational velocities (Table [5|) , compared to single field G- 
type stars, we have investigated such models. 

Models with reduced mixing-length parameters, includ- 
ing non-gray atmospheres and an improved equation of 
state, have been published by Baraffe et al. (jT998j3, but 
only for solar metallicity (up to 1.4 Mq) and [Fe/H] = —0.5 
(up to 1.0 Mq). 

Guided by a comparison between V636Cen and these 
models, we have calculated dedicated models for the ob- 
served masses and [Fe/H], tuning the mixing length param- 
eter individually for the components. We have adopted the 
code by Claret PIMIIMISI [MMl [MI7|) . which assumes an 
enrichment law of y = 0.24 ± 2.3Z together with the solar 
mixture by Grevesse & Sauval (11998) , leading to {X,Y,Z)q 
= (0.704,0.279,0.017). The observed [Fe/H] = -0.20 then 
corresponds to {X,Y,Z) = (0.724,0.264,0.012). We note 
that the envelope mixing length parameter needed to re- 
produce the Sun is I /Hp = 1.68. 

Four different cases, listed in Table [HI will be discussed 
here. Fig.[8]shows temperature-radius comparisons between 
V636 Cen and models with different mixing length parame- 
ters, calculated for the component masses. For the observed 
[Fe/H], the temperature of the primary component is fitted 
well by the I /Hp = 1.4 model, whereas 1.0 is adequate for 
the secondary (case 1). As seen from the age-radius diagram 
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Fig. 8. V636 Cen compared to Claret models for the ob- 
served masses, calculated with different mixing length pa- 
rameters. The four cases listed in Table [5] are shown. The 
circles on the evolutionary tracks mark model ages of 1.0, 
2.0, and 3.0 Gyr. Case 1 = full drawn, thick; Case 2 = 
dashed; Case 3 = full drawn, thin, open circles; Case 4 = 
dotted. 

in Fig. [9l these models also predict practically identical ages 
for the components. Changing the mixing length param- 
eters by just -1-0.1 (case 2), i.e. considering models which 
still fit the observed masses and temperatures within errors 
(Fig. [5]), lead to significantly higher ages, which are, how- 
ever, no longer identical for the components. To illustrate 
the effect of the abundance uncertainty ±0.08 dex, models 
for [Fe/H] = —0.12 have been included in the two figures 
(case 3). For this abundance, mixing length parameters of 
about 1.6 (primary) and 1.1 (secondary) are needed to fit 
the observed masses and temperatures. Ages determined 
from the masses and radii become even higher, and they 
differ for the components. Finally, Claret models adopting 
a solar I /Hp are, like the Y'^ models, unable to fit V636 Cen 
(case 4). 

Table [5] summarizes the ages derived from the different 
model cases. Effects from the mass uncertainties are negli- 
gible and have not been included; see Fig. O In conclusion, 
only the case 1 models reproduce all observed properties of 
the components of V636Cen well at a common age (1.35 
Gyr). This implies that I /Hp can be tuned fairly precisely, 
provided accurate dimensions and abundances are avail- 
able. We finally note, that according to the Claret models 
both components have pronounced convective envelopes, 
starting at 0.78 Rp and 0.73 Rs- 

7.3. Comparison with other binaries 

Currently, nine other main sequence eclipsing binaries with 
at least one component in the 0.80-1.10 Mq mass interval 
have well-established (« 2% or better) masses and radiPIF°l. 

® see e.g. Ihttp : //mmm . astro . keel e . ac . uk/"^ j kt7| 

Two additional binaries in this meiss interval, the pre-main- 
sequence system V1174 Ori (Stassun et al. l2004p . and V432 Aur, 




0.80 - 



0.75 r","' , , , I , , , I , I , , , , I , , , , I , , , , I , , , , I , , , , ' 
0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 
Age (Gyr) 

Fig. 9. V636Cen compared to Claret models. The four 
cases listed in Table [HI are shown. The curves illustrate 
model radii as function of age for the component masses 
(upper: primary; lower: secondary). Case 1 = full drawn, 
thick; Case 2 = dashed; Case 3 = full drawn, thin; Case 4 
= dotted. The horizontal full drawn lines mark the observed 
radii of the components with errors (dotted lines). The ver- 
tical line marks the common age (1.35 Gyr) predicted by 
the case 1 models. 

Furthermore, the A and B components of of a Cen can be 
added, since their masses, radii, temperatures, and abun- 
dances are known to high accuracy (e.g. Eggenberger et al. 

These 1 1 binaries are listed in Table [TO] and shown in 
Fig. [TO] together with solar composition isochrones. We 
note, that V636 Cen and a Cen are the only systems in this 
sample with [Fe/H] measured spectroscopically. With one 
exception, VZ Cep, the systems fall in two groups: 

1. ForCGCyg, CVBoo, ZZUMa, FLLyr, and V1061Cyg, 
like for V636 Cen, the Y^ models predict higher ages for 
the secondary components than for the primary compo- 
nents, resulting from larger radii than predicted. Also, 
the effective temperatures of the secondary components 
are lower than predicted. Their orbital periods are short 
(0'?6-2^8), and the rotational velocities of their compo- 
nents are high (20-73 kms~^). Intrinsic variations are 
present, and some of the binaries are known, strong X- 
ray sources. Information on Ca ii emission is unfortu- 
nately scarce. A few further systems with similar char- 
acteristics, but less accurate dimensions (RT And, BH 
Vir, UV Leo, and UV Psc), could be added to this group 
(Clausen et al. [T999)) . 

2. For HSAur, RWLac, a Cen, and V568Lyr (V20 in 
NGC6791), the observed masses, radii, and tempera- 
tures agree much better with those predicted by the 
solar-scaled Y^ models. They have fairly long orbital pe- 
riods (9'?9-79.9 yr), and the rotational velocities of the 
components are below 6 kms^^. Popper et al. ( 11986P 
saw no Ca ii emission but found evidence for intrinsic 

which has an evolved component (Siviero et al. I2004|) . have not 
been been included in the comparison. 
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Table 10. Properties of well-studied main sequence binaries (e.g. footnote 9) with at least one component in the 
0.80-1.10 Mq mass interval. 





PGriod 


M IMr^ 




lug Lij UQ 


v sin 


V CXiL i-CL Vji.1.1 \j y 


TT Inpf T,v^ 












(kms~^) 
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" If measured rotational velocities are not available, synchronous values (s) are given 
' X-ray luminosities are based on Voges et al. (|1999l I2000|l and Motch et al. (|1997|l 



variability of HSAur and/or the KO comparison star 
used (HD49987). The other three binaries are reported 
to be stable; no information on Ca 11 emission is avail- 
able. None of the four is a known strong X-ray source. 
According to models, RWLac is somewhat metal defi- 
cient, and its age is about 11 Gyr (Lacy et at. I2005p . 
V568 Lyr has an age of about 8 Gyr and its metal con- 
tent is about 2.5 times solar (Grundahl et al. I2008p . 



For the exception, VZCep, the observed masses and 
radii are well reproduced by the models, but the ob- 
served temperature of the secondary component is about 
250 K lower than predicted (Torres et al. I2009p . Therefore, 
it is also less luminous than predicted, which is in general 
not the case for the rest of the sample. VZ Cep is a strong 
X-ray source, and as for the binaries in group 1, the or- 
bital period is short and the components rotate fast. On 
the other hand, no intrinsic variations are reported. 

From a comparison of V1061Cyg and FLLyr (group 
1) against RWLac and HSAur (group 2), Torres et al. 
(|2006p identified chromospheric activity as the likely cause 
of the radius and temperature discrepancies. Our sample 
strengthens this suggestion, not at least through a compar- 
ison between V636Cen and V568Lyr, whose components 
have nearly identical masses. 



7.4. Mixing-length calibrations 

It is worth linking our sample to the ongoing discussion on 
the mixing-length parameter for stars in the 0.7-1.1 Mq 
range, where stars have convective envelopes but are not 
fully convective. From 2D radiation hydrodynamics (RHD) 
calculations, Ludwig et al. (1999) calibrated the mixing- 
length parameter for solar-type stars (Toft 4300-7100 K; 
logg 2.54-4.74) and found a moderate but significant vari- 
ation of I /Hp. On the main sequence, a decrease with in- 
creasing temperature, and thereby mass, from about 1.8 to 
about 1.3 was obtained. On the other hand, Yildiz et al. 
(|2006p and Yildiz (2008,) found a significant increase with 
mass from model fits to four Hyades binaries, FL Lyr, and 
V442 Cyg. As an illustration, I /Hp values calculated from 
the fitting functions by these authors are listed in Table [TTl 
no scaling, e.g. to the Sun, has been done. Finally, including 
seismic constraints, several studies of a Cen, most recently 
by Eggenberger et al. pM|) and Yildiz (PDH7|) . find that 
I /Hp is higher for the less massive B component than for 
the A component. Eggenberger et al., who use the seismo- 
logical data by Carrier & Bourban (|2003p . obtain 1.83 (A) 
and 1.97 (B), and Yildiz, who include the data by Kjeldsen 
et al. (120051) . get 1.64 (A) and 2.10 (B). 

The fact that our sample can be divided in two different 
groups seems to suggest that the behaviour of the mixing 
length parameter may be governed by two different effects: 

— a slight decrease of I /Hp with increasing tempera- 
ture/mass for inactive main sequence stars, as predicted 
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Fig. 10. Main-sequence binaries with one or both compo- 
nents in the 0.80-1.10 Mq interval liaving well-established 
masses and radii (2% or better); see Tabel [TUl The compo- 
nents of V636 Cen are connected by a solid line. The dottes 
lines are isochrones (0.5 Gyr and 1-10 Gyr, step 1 Gyr) 
for Z — 0.0181, i.e. solar composition (Demarque et al. 
[2004]) . 

Table 11. Mixing- length parameters {1/ Hp) for the binary 
components (Table [TU)) calculated from the temperature- 
gravity fitting function by Ludwig et al. (1999)) and the 
mass fitting formula by Yildiz et al. (|2006p . 
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1.58 
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1.60 
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1.58 


1.61 


2.06 


1.78 
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1.59 
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1.42 
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VZ Cep 


1.39 


1.59 


2.34 


2.07 



" For comparison, the values calculated for the Sun are in- 
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by Ludwig et al. (|1999p and found for a Cen and possi- 
bly also V568Lyr (Grundahl et aL [M]5|) . and 
— a decrease of I /Hp for active stars compared to inactive 
ones of similar mass, as seen for group 2. 

If this is the case, the result by Yildiz et al. (|2006p and 
Yildiz (I2008P might then be questioned, since at least two of 
the binaries in their samples, the Hyades binary VSlSTau 
and FLLyr, have active component (s). 

It would be of interest to extend the systematic study 
of the evolution of low-mass stars and brown dwarf eclips- 
ing binaries by Chabrier et al. (|2007p to higher masses. 
For masses up to 0.8 Mq, they found that the effects of 
decreasing I /Hp, i.e. reducing the convective efficiency, in- 
creases with mass, leading to lower effective temperatures, 
luminosities, central temperatures, and nuclear energy pro- 
duction. In parallel, magnetic spots covering a significant 
fraction of the stellar surfaces cause similar changes. They 
adopted a phenomenological approach, and more complete 
models, which directly couple rotation, convection, and dy- 
namo magnetic fields are also desirable; see e.g. D'Antona 
et al. pOOOp for pre-main-sequence models. 

8. Summary and conclusions 

From state-of-the-art observations and analyses, precise 
(0.5%) absolute dimensions have been obtained for the to- 
tally eclipsing solar-type system V636 Cen. A detailed spec- 
troscopic analysis yields a Fe abundance relative to the Sun 
of [Fe/H] = —0.20 ± 0.08 and similar relative abundances 
for Si, Ca, Ti, V, Cr, Co, and Ni. 

The 0.85 Mq secondary component is found to be mod- 
erately active with starspots and Ca ii H and K emission, 
and the 1.05 Mq primary shows signs of activity as well, 
but at a much lower level. Apsidal motion (C/ = 5 270 ±335 
yr) with a 40% relativistic contribution has been detected 
for the eccentric orbit (e = 0.135 ± 0.001), and the inferred 
mean central density concentration coefhcient, log(A:2) = 
— 1.61 ± 0.05, agrees marginally with model predictions. 
The measured rotational velocities, 13.0 ± 0.2 (primary) 
and 11.2±0.5 (secondary) kms~^, are in remarkable agree- 
ment with the theoretically predicted pseudo-synchronous 
velocities, but are about 15% lower than the periastron val- 
ues. 

We have shown that stellar models with solar-scaled 
mixing length parameters are unable to match the com- 
ponents of V636Cen at identical ages. At the age of the 
primary, the secondary component is ~ 10% larger than 
predicted. Also, the components are ~ 200 K (primary) 
and ^ 400 K (secondary) cooler than predicted. 

However, models adopting significantly lower mixing- 
length parameters can remove these discrepancies. For the 
observed [Fe/H], Claret models for l/Hp — 1.4 (primary) 
and 1.0 (secondary) reproduce the properties of the com- 
ponents of V636 Cen well for a common age of 1.35 Gyr. 

Currently, 10 other solar- type binaries have well- 
estqblished dimensions, although spectroscopic abundance 
results are lacking for most of them; see Fig. [TO] and 
Table [TO] With one exception, VZ Cep, they fall in two 
groups: Four long period, slowly rotating inactive systems, 
which seem to be fitted reasonably well by solar-scaled 
models. And five systems, which exhibit intrinsic varia- 
tion, and have short orbital periods and high rotational 
velocities. Like V636 Cen, they can not be reproduced by 
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solar-scaled models. Therefore, the sample strengthens the 
suggestion by Torres et al. (|2006p . based on fewer systems, 
that chromospheric activity, and the associated changes of 
the convective envelope, is the likely cause of the radius and 
temperature discrepancies. The strongest case is perhaps a 
comparison between V636 Cen (active) and the newly stud- 
ied V568Lyr in NGC6791 (inactive), whose components 
have nearly identical masses. 

Finally, the clear division of our sample in two groups 
might suggest that two different effects influence the mixing 
length parameter for 0.7-1.1 Mq main sequence stars: a) a 
slight decrease of I /Hp with increasing temperature/mass 
for inactive main sequence stars, as found by e.g. Ludwig 
et al. (|1999|) and seen for a Cen and possibly also V568 Lyr, 
and b) a decrease of I / Hp for active stars compared to in- 
active ones of similar mass. 

More well-studied eclipsing binaries, as well as more 
complete stellar models, are needed to fully understand and 
describe envelope convection in terms of effective temper- 
ature, gravity, mass, age, rotation, dynamo magnetic field, 
activity level etc. We are presently studying several systems 
with main sequence components in the 0.5-1.2 Mq range, 
including V1123Tau, V963Cen, ALDor, QRHya, KYHya, 
ALAri , UW LMi, EWOri, and NYHya; see also Clausen 
et al. (|200ip . They exhibit various levels of activity, from 
none to modest, and should therefore provide valuable new 
insight into this matter. Spectroscopic abundance determi- 
nations will be part of the analyses. 
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Table A.l. Radial velocities of V636 Cen and residuals from the final spectroscopic orbit (in units of kms ^). a is the 
internal velocity error. 
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